Derivation of non-Markoffian transport equations for trapped cold atoms in 

nonequilibrium thermal field theory 
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The non-Markoffian transport equations for the systems of cold Bose atoms confined by a exter- 
nal potential both without and with a Bose-Einstein condensate are derived in the framework of 
nonequilibrium thermal filed theory (Thermo Field Dynamics). Our key elements are an explicit 
particle representation and a self-consistent renormalization condition which are essential in thermal 
field theory. The non-Markoffian transport equation for the non-condensed system, derived at the 
two-loop level, is reduced in the Markoffian limit to the ordinary quantum Boltzmann equation 
derived in the other methods. For the condensed system, we derive a new transport equation with 
an additional collision term which becomes important in the Landau instability. 

PACS numbers: 03.75.Kk, 05.20.Dd, 05.30.Jp 



I. INTRODUCTION 

The systems of trapped cold atoms are ideal for study- 
ing quantum many-body theories such as quantum field 
theory and thermal field theory. They are dilute and 
weak-interacting, so theoretical calculations can be com- 
pared with experimental results directly. 

Since the realization of Bose-Einstein condensates 
0) 0) 9 > the formation and grow of condensate [1] , the 
thermal shift of the energy spectrum [B[ , and many other 
intriguing phenomena have been observed with good ac- 
curacy, and offer opportunities to test quantum many- 
body theories in both equilibrium and nonequilibrium. 

In the aim of describing the kinetics of the trapped cold 
atom system, a number of theoretical approaches have 
been proposed such as the methods of the quantum Boltz- 
mann master equation [|| 0], the quantum Boltzmann 
equation with the local density approximation [8|, |9| , the 
closed path (CTP) formalism [l(J, HH, an d the effective 
Hamiltonian method in Thermo Field Dynamics (TFD) 
12fl. These are in good agreement with the experiments 
2, |J] ■ They however are based on a phase-space distri- 
bution function, and the energy spectrum is not quan- 
tized. This implies that the discussions of the particle 
representation or the diagonalization of the Hamiltonian 
are absent, while they are essential for the quantum field 
theory. 

There are two nonequilibrium extension of the ther- 
mal field theory, i.e., the closed time path formalism 



'Electronic address: nakamura@aoni.wascda.jp 
i Electronic address: 



tomoka@fuji.waseda.jp 
mine@waseda.jp 



+ Electronic address: 
^Electronic address: okumura.masahiko@jaea.go.jp 
^Electronic address: yamanaka@waseda.jp 



[13] and TFD 14]. The CTP formalism is widely used. 
But we employ the TFD formalism in this paper, be- 
cause the concept of quasi-particle picture is clear even 
in nonequilibrium situations there. In TFD, which is a 
real-time canonical formalism of quantum field theory, 
thermal fluctuation is introduced through doubling the 
degrees of the freedom, and the mixed state expectation 
is replaced by an average of a pure state vacuum, called 
the thermal vacuum. 

It is crucial in our formulation of TFD to construct the 
interaction picture. In quantum field theory, the choice 
of unperturbed Hamiltonian and fields is that of quasi 
particle picture, and concrete calculations are possible 
only when a particular unperturbed representation, or 
a particular particle picture, is specified. One does not 
know an exact unperturbed representation beforehand. 
Taking plausible representations, parameterized by some 
parameters, we calculate the propagators of the Heisen- 
berg fields and require some conditions on them, called 
self-consistent renormalization conditions, which pick up 
a self-consistent representation and determine the param- 
eters. The renormalized mass and coupling constant are 
such examples in quantum field theory. We construct 
a quasi particle picture in the doubled Fock space in 
TFD, defining quasi particle operators which diagonal- 
ize the unperturbed TFD Hamiltonian. In nonequilib- 
rium case a time-dependent number distribution is in- 
troduced as a unknown parameter, and a self-consistent 
renormalization condition derives a n eq uation for it, i.e., 
the quantum Boltzmann equation [141 ]. Moreover, the 
non-Markoffian extension of the self-consistent renormal- 
ization condition is also proposed [151 ]. However, the ex- 
tension and application to intrinsically inhomogeneous 
systems and condensed ones have not been established. 

In this paper, we derive the non-Markoffian quantum 
transport equations for cold atoms in a confining po- 
tential both without and with a condensate from the 
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nonequilibrium TFD formalism [Til Il5| . We confirm 
that our non-Markoffian transport equation for the non- 
condensed system is reduced to the ordinary quantum 
Boltzmann equation derived in the other methods when 
the Markov approximation is applied. For the condensed 
system, we find that the non-Markoffian equation con- 
tains an additional collision term which is overlooked in 
the other methods. This term vanishes in the equilib- 
rium limit if there is no Landau instability, but remains 
non- vanishing to prevent the system from equilibrating if 
there is Landau instability. Thus our transport equation 
with the additional term (we call it the triple produc- 
tion term) and the other ones without it predict defi- 
nitely different behaviors of the unstable system. This 
difference is traced back to different quasi particle pic- 
tures in the respective theories. Although we only con- 
sider in this paper the systems of trapped Bose atoms, 
our formulation of nonequilibrium TFD can be extended 
straightforwardly to the trapped systems of Fermi or 
multi-component atoms. 

This paper is organized as follows. We briefly review 
the formulation of nonequilibrium TFD in Sec. II. In 
Sec. Ill and IV, the non-condensed and condensed sys- 
tems are considered, respectively. We formulate each in- 
teraction picture corresponding each quasi particle pic- 
ture, diagonalizing the free (unperturbed) Hamiltonians 
of TFD. The tensor form [16j |. which makes the diagram- 
matic calculation very simple, is introduced for the non- 
condensed system, and is extended to the condensed sys- 
tem. Applying the self-consistent renormalization condi- 
tion proposed by Chu and Umezawa [Hj], we construct a 
systematical method to obtain the transport equation for 
the trapped systems. Section V is devoted to summary 
and discussions. 



II. NONEQUILIBRIUM TFD FORMULATION 

Here we briefly review the formulation of nonequilib- 
rium TFD [3. 

In TFD, every operator A gets its tilde conjugation 
pair A , which is related to the ordinary (non-tilde) op- 
erator by the following tilde conjugation rules: 



(AB)~ 


= AB, 


(1) 


c 2 B)~ 


= c\A + c*B , 


(2) 




= A, 


(3) 


(At)~ 


= A\ 


(4) 


|0}~ 


= 10), 


(5) 


(Of 


= (0|, 


(6) 



where c\ and ci are arbitrary c- numbers, and |0) and 
(0| are the thermal vacua. The Hamiltonian of TFD, 
which should generate the time translations of both non- 
tilde and tilde operators, is not the ordinary Hamiltonian 
H but the hat Hamiltonian H = H — H . The time 



independence of the thermal vacua requires the minus 
sign in front of H . 

The construction of the interaction picture is crucial, 
because the choice of an interaction picture corresponds 
to that of a quasi particle picture. Suppose the bosonic 
a^-operators in the interaction picture, representing the 
renormalized quasi particle with a quantum number t. 
They are related to the ^-operators (called the repre- 
sentation particle operators), which annihilate the time 
independent thermal vacuum £^|0) = 0, through the ther- 
mal Bogoliubov transformations 

a%{t)=B^»(t)$(t), (7) 

a^(t)=e;(t)Br(t), (8) 

e t {t)=B>?{i)am, (9) 

e i (t) = a?(t)B7 1 '» v (t). (10) 

Here we introduce the thermal doublet notations 

a^t^Y , a"=( t -~a) V , (11) 

=(!)"' ? = @ -eT. (12) 

and the thermal Bogoliubov matrix 

sr(t) = A + »,(() -»,(t)Y" (13) 

It is important to take the above particular form of the 
thermal Bogoliubov matrix, as one calls a = 1 represen- 
tation [7JJ , which enables us to make use of the Feynman 
diagram method in nonequilibrium systems 17]. The 
number distribution ni(t) is given by 

n t (t) = (0\4(t)a e {t)\0) , (15) 

and its time dependence is determined later. 

The unperturbed Hamiltonian for the ^-operators 
should be diagonal, consistently with the time indepen- 
dence of the thermal vacuum. So the time dependence 
of ^-operator in the interaction picture should be in 
the form $(t) = ^e-^'^^W , generated by the free 
Hamiltonian Ho(t) — J2 £ wi(t)£%{t)£%(t) . Throughout 
this paper h is set to be unity. Note that cue generally 
depends on time because of the time dependent energy 
renormalization. In this paper, we take time independent 
u>£, assuming that the energy shift is negligible in the 
leading order of perturbation, i.e., £^(i) = ^e^' luJet and 
Hq = ^2^uji^(t)^{t). The unperturbed Hamiltonian 
for the a^-operators is not Hq, but Hq(t) — Hq — Q(t) 
with the thermal counter term Q(t) 

Q(t) = i^ht(t)aZ(t)Q :fpa£(i) (16) 
= -<£>(*)#(*) (S jV^W, (17) 
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caused by the i-dependence of ni{t). 

The field operator ip(x) is expanded with a complete 
set {ue(x)} as 



(18) 



where x = (x,t). The unperturbed and full propagators 
for ip and £ are defined by 

A^(x,x') = -i(0\T[^(x)^ u (x')}\0), (19) 

G^(x,x') = -i(0\T[^(x)^(x')]\0) , (20) 

d%,(t,t') = -i{0\T[$(t)&(t')]\0), (21) 

g^f) = -i(0|T[^(t)^,(t')]|0) , (22) 

respectively, which are related to each other as 



A" v (x,x') 



G^(x,x') 



'u e (x)B_ 



xd'tr {t,t')By{t')uUx'), (23) 



xgHf {t,t')BV{t')ul,{x'). (24) 



The subscript h denotes a quantity in the Heisenberg pic- 
ture. While the unperturbed propagator d has a diagonal 
structure 



d%,(t,t') = 5, 



-i6{t-t') 
i6{t' - 1) 



-ioji{t—t ) 



(25) 

the full propagator g has an upper triangular structure 
in general, that is, gj}(t,t') ^ and gjj,(t,t') = [lij]. 
This is because ^ and fjj identically annihilate the bra- 
vacuum in the a — 1 representation while £h and £h 
do not generally annihilate the ket-vacuum. It is also 
shown that g 11 and <? 22 are a retarded and an advanced 
functions, respectively, as d 11 and d 22 are. 

According to the Feynman method, one can calculate 
the full propagator with the interaction Hamiltonian in 
the interaction picture, 



Ht — H u 



Q. 



(26) 



with fljnt = H — Hq . Possible renormalization counter 
terms are suppressed below for simplicity. 



III. TRAPPED BOSE ATOMS IN 
NON-CONDENSED SYSTEM 

In this section, we derive the transport equation for 
the system of cold Bose atoms without condensate in 
nonequilibrium TFD. 

We start with the following Hamiltonian to describe 
the trapped dilute Bose atoms, 



H= d 6 



(27) 



where m, V(x), n, and g represent the mass of an atom, 
the trap potential, the chemical potential, and the cou- 
pling constant, respectively. The bosonic field operator 
4>{x) obeys the canonical commutation relations 

W{x)^{x')]\ t=tl =8{x-x'), (28) 
[4>{x),i>(x')]\ t=t , = [^(x), VV)]| t=t , =0. (29) 

We expand the field operator ip(x) as in Eq. (fT8|) . us- 
ing the solutions of the following eigenequations, {ug(x)} 
with the eigenvalues {^e}, 



(30) 



— V 2 + V(x) - /i ) u e (x) = w t ut{x) 



The annihilation- and creation-operators ae and a\ diag- 
onalize the free Hamiltonian part Ha 



at ■ 



H = Jd 3 x^ ("^V 2 +V(x)- fiji> = J2"t4 

(31) 

We apply the formulation of nonequilibrium TFD in 
the previous section to the present system: Each degree 
of freedom is doubled, the time dependent thermal Bo- 
goliubov transformation is introduced in the interaction 
picture, and the total Hamiltonian H is divided into the 
unperturbed and interaction parts, Hq and Hi . Then 
the full propagator is calculated in the Feynman diagram 
method. 

The self-consistent renormalization condition on the 
full propagator thus obtained, which is extended from 
the self-consistent on-shell renormalization condition in 
the ordinary quantum field theory, is already proposed 
[lldllli as 



<?«(M)=0. 



(32) 



It provides the transport equation which determines the 
temporal evolution of the unperturbed number distribu- 
tion ni(t) . Following the Dyson equations G — A + AEG 
or g = d + dSg, we obtain 

9U> (*> t ')=Y. j dSClS ' »i™(*> S )^«>' ( S ' S '^fn'i' ( S '> O » 
mm' 

with the self-energy 

Z^(x,x') =Y,ni{x)BJ x '^' \t) 

x S&? (t, t')B$ u (t')u* e , (x') , (34) 

anv/. +/\ _ fS}i,(t,t') S} 2 ,(t,t')\ , . 

b u \t,t)-y Q S } 2 (t,t')) ■ {6b > 

To illustrate how the transport equation follows from 
the renormalization condition, we approximate the full 
propagators g 11 ' 22 ^ in Eq. (|3"3"|) by the unperturbed ones 
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d 11 ( 22 ) and divide the self-energy S into a loop contri- 
bution Si oc ,p and a contribution of the thermal counter 
term Sq, 



S^ u ,(t,t') = -ih e (t)5wS(t-t r ) 



Then we have 



1 





(36) 



9u (*) t) = -i ds 



2Rc 



fli(s) 

ds e luJe<,s s ^ Sli i oop (s, s ) 



(37) 



and the renormalization condition 
lowing transport equation 



implies the fol- 



n e (t) = 2Re / dse^-^S^Jt^s) 



(38) 



In view of this, we are going to calculate the self-energy 
perturbatively to obtain the transport equation in the 
leading order. Before that, we introduce the following 
tensor form [Til ] which makes the representations and 
calculations of propagators and self-energies much more 
concise 



si (if ^ = 1) 
s 2 (if /i = 2) 



(39) 



For instance, the following matrix appearing in the un- 
perturbed propagator is expressed in the tensor form as 



B(ty 



1 





B(f) 



l + n(f) -n(i') xA "' 
l + n(f) -n{t'\ 

" r l + n{t'Y 
-n(t') 



(40) 
(41) 



Thus, the unperturbed propagators can be written in the 
tensor form as 



-id we 



-iujc(t' — i) 



d%,(t,t>) 



(42) 



m{t) 

+ n e (t) 



-iujg (t—t ) 



+ o(t' -t)\ 1 



(43) 



with the sign factor, e 1 — 1 and e 2 = — 1 . 

As an example of manipulating products of the unper- 
turbed propagators in the Feynman diagram calculation, 




FIG. 1: Two- loop self-energy diagram. 



we give the following manipulation, 



' (J || ) 



Beit)' 1 (J °W*0 



l + n e (t')V (lV (l + ne(t')V 



-n t {t') 



-n v {t') 



11" Ul + ne(t'))(l + n e >(t')) 
if \ nt(f)n v (?) 



(44) 
(45) 



where the indices \i and v are not summed over. This 
kind of manipulation will be used below in the calcula- 
tions of the self-energies. 

We focus on the two-loop self-energy indicated in 
Fig. [U the leading loop diagram which makes g 12 (t,t') 
nonzero, 

V&vM = -2g 2 e^"A^(x,x')A^(x,x')A^(x',x) . 

(46) 

The sign factors arise because of the definition of V> and of 
the particular form of the interaction Hamiltonian Hi = 
Hi — Hi . Using the tensor form, we can rewrite the 
self-energy as 

T,£ op (x,x') = -2ig 2 e v ^ e -*(«« l +Wi a -«i,)(t-0 

x u ll (x)u£ 2 {x)u* l3 (x)u* ti (x')u* t2 (x')ui s (x') 

x 9(t -*){}}" + + 

Mil 1 n ei rii 2 (l + n i3 ) 



\(1 + nfj(l + ne 2 )n is 



min(t.t') 

(47) 



where the subscript min(i, t') denotes the time arguments 
of n. So that 



where 



n^m^l + nt 3 ){\ + n e ) 

- (1 + n^)(l + n e2 )n i3 n e 



min(t,t') 



d x {x)ui 2 {x)u*« {x)u} (x) 



(48) 
(49) 
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Thus the transport equation at the two-loop level is de- 
rived 

ht(t) = 4g 2 Ref ds e-^+^-^-^X*-*) C llW 

x n^nizil + ne 3 )(l + n e ) - (I + n ei )(l + ni 2 )ni 3 n e . 

(50) 

The transport equation of Markoffian type correspond- 
ing to Eq. (l50| has already been derived by Chu and 
Umezawa [15[ for the homogeneous system. We have 
here extended their result to the inhomogeneous system 
due to the confining potential. 

The exponential term e - '^!" 1 "^ - ^ in 
Eq. (j5"0)) signifies the energy conservation. To see that 
explicitly, we perform the time integral in Eq. (I50|) as 



ht{t) = 4<7 2 Re 



J — OO nun 



J(ut 1 +U£ 2 -<+>e 3 —Ul— 



x Ci 1 e 2 -£ 3 £ Ri 1 £ 2; e 3 e(t) (51) 
(52) 



where 



Rl x l 2 ;i 3 l = n<,n< 2 (l + »f s )(l + nt) 

- (1 + n^)(l + ni 2 )ni 3 n e , (53) 



dR e 



R, 



dt 



(54) 



ft is seen from this Lorentzian form that the energy is 
conserved with the width T which is small when the tem- 
poral change of n is slow. 

To confirm the correspondence of our transport equa- 
tion with those derived in the other methods, we approx- 
imately replace ng i (s) in the right hand side of Eq. ([50]) 
with n^(t), in other words, take the limit r — * or the 
Markoffian limit, with the assumption that the system is 
close to the equilibrium and the evolution is sufficiently 
slow, 



nt 



(t) = Aug 2 2^ 8{wi 1 + uj£ 2 - 0Jt 3 - u>t) Ct x t 2 ;i 3 i 



x n tl ni 2 (l + 71^,) (1 + nt) - (1 + n tl )(l + ni 2 )n f3 n e 

(55) 

The equation has a form of the ordinary quantum Boltz- 
mann equation and is consistent with the one which has 
been obtained from the on-shcll renormalization condi- 
tion for the homogeneous system in nonequilibrium TFD 
[l4j |. However, the naive approximation is not valid for 
a trapped system with the discrete energy spectrum u>t. 
The energy conservation forced by the delta function is 



too strict to allow any energy exchange of particles and, 
consequently, any time evolution. One way to avoid this 
difficulty is to apply the coarse graining treatment [20j . 
For instance, in the case of the harmonic trap potential, 
S(w tl +u>i 2 -u>t 3 -uJt) was replaced by 8 Uil+U)l3tUts+ut /to 
with the trap frequency f2 [2l[ , 

47T<? 2 



x n ei n£ 2 (l + n i3 )(l + nt) - (1 + n fl )(l + ne 2 )n£ 3 n e 



(56) 

Note that such simple replacement is valid only for the 
harmonic trap system whose energy-level spacing is uni- 
form. 

Another way is a semi-classical formulation using the 
phase-space distribution function n(x,p, t). Then the 
particle energy is no longer the discrete one but the local 
continuous one u>(x,p) = p 2 /2m + V(x) — (i, for which 
the difficulty mentioned above does not arise. Although 
the transport equations, derived in this manner by sev- 
eral authors 0, 0, EH , are in good agreements with the 
experiments of evaporative cooling and formation of con- 
densate, they are classical and can not describe quantum 
fluctuations fully, and the particle picture is not explicit. 
In the next section we will show the case where the in- 
stability of condensate exists and the quasi particle spec- 
trum plays a crucial role, for which the semi-classical 
treatment is not valid. 



IV. TRAPPED BOSE ATOMS IN CONDENSED 
SYSTEM 

In this section, we consider the situation in which a 
condensate exists and derive the transport equation. The 
field operator ip(x) is divided into a classical part £(cc) 
and a quantum part <p{x), reflecting the existence of the 
condensate. In a fully nonequilibrium situation, the or- 
der parameter £ = (0 1 -0 1 0) should be time-dependent, 
however we consider only a situation near the equilib- 
rium and assume the time-independent order parameter 
throughout this paper. 

The Hamiltonian l|27p is written as 



H — H + H int , 



(57) 



where 
Ho 

H\„t 



d 6 x 



2m 



+ V{x)-n + 2g\Q{ X )\ 2 \y 



+ |(c* 2 ( a # 2 + cW< 2 ) 



, (58) 



d A x 



(59) 
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and the first order term of tp(x) vanishes since C( x ) is re_ 
quired to satisfy the following Gross-Pitaevskii equation 
22] at the tree level 



2m 



+ V{x)-n + g\{(x)\ J ){{x) = 



(60) 



Next, we briefly review the Bogoliubov-de Gennes 
(BdG) method which diagonalizes the free Hamiltonian 
H . The BdG equations are simultaneous eigenvalue 
equations given by [H [H HH] 



Tyi(x) = tx>zyi{x) . 
Here the doublet notation is introduced as 



Vl 



T = 



L M 

-M* — L 



where 



L 

M=gC(x) 



2m 

'2. 



+ V(x)-n + 2g\{(x)\ 



(61) 

(62) 
(63) 



(64) 
(65) 



It is known that the BdG equations have the zero eigen- 
value mode whose treatment needs attention. We disre- 
gard the zero mode for simplicity though it can be in- 
cluded consistently [12, HH, because the scope of this 
paper is confined to the time-independent order param- 
eter and the dynamical effects of the zero mode is not 
dominant then. In addition, the eigenvalues can be com- 
plex since the operator T is non-Hcrmitian. The condi- 
tion for the emergence of complex eigenvalues in the BdG 
equations has been studied both numerically [13, [H, [29j] 
and analytically [3(], HH, HH , and the quantum field the- 
oretical formulation has also been discussed [341 ]. The 
emergence of complex eigenvalues implies the dynamical 
instability of the system, and a drastic temporal change 
of the order parameter occurs then, which is out of our 
present formulation. In this paper, we consider only the 
case where no complex eigenvalue emerges. 

Eigenfunctions belonging to the non-zero real eigenval- 
ues can be orthonormalizcd under the indefinite metric 



as 



Jd 3 x yl(x)a 3 y e >{x) = 5 U 
Jd 3 x z\(x)o 3 ze'(x) = —I 
Jd 3 x y\{x)a 3 z t >{x) = , 



(66) 
(67) 
(68) 



with i-th Pauli matrix cr.;. The function ze, defined by 
ze — <J\y\, is an eigenfunction belonging to —uit, if ye 
is an eigenfunction belonging to tog. It is convenient to 



rewrite the orthonormal conditions (|66)) - (|68|) with the 
2x2 matrix form as 

Jd 3 x W t (x) W^(x) = 5a> , (69) 

where 

W t (x)=a 3 ( V \[ X Aa a , 
\z}(x)J 

Wf\x) = (y e {x) z e (x)) . 
The completeness condition, 

[vt{x)y\{x') - z t {x)z\{x')\ = a 3 S(x - x') , (72) 



(70) 
(71) 



can be expressed as 

^2 W^ix) Wi(x') = 8{x - x') . 



(73) 



and then the field operators are expanded in the doublet 
form as 



V a {x)=Y,Wi 1 ^{x)b^t). 



>(x)=Y,b?(t)Wf(x), 



where 



bf = {b\ -b t y 



(74) 
(75) 

(76) 
(77) 



The operators bp, satisfy the canonical commutation rela- 
tion [be , b\,] = Sw, and diagonalizes the free Hamiltonian 

(El 



H Q = l - jd 3 X <p a 7 r - ; (.r) 



Y^"ib\l 



(78) 
(79) 



The operators be annihilate the Bose-Einstein con- 
densed vacuum and the operation of the creation opera- 
tors be on the vacuum constructs the Fock space at zero 
temperature. Therefore, to treat this system in nonequi- 
librium TFD, we double the degrees of freedom as follows 



(bl -b 



b: 



= fe f 



if 

6) Be 



(80) 
(81) 



where the thermal Bogoliubov matrix is defined in 
Eqs. (fl"3|) and (TT4")) with the quasi particle distribution 
n e(t) — (0\be(t)be(t)\0) . Note that the operators who an- 
nihilate the thermal vacuum are not the 6-operators but 
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the ^-operators. The combination of the two transforma- 
tions, £ into b and b into (p, involves the 4x4 transfor- 
mations, 



fb t \ 

b\ 
b\ 

w 

M 
W 



1 ne 


\ 






1 l+ni 










1 + n e 1 








7l£ 1 / 







£ 





4 \ 


(h\ 


vi 


4 


b\ 




4 




V v 2 i 


4 J 


W 



(82) 



(83) 



(84) 



where the blank elements denote zero. It is convenient 
to introduce the quartet notations for bg as follows 

b\ 

\btj 

(&? fM")^^ -fc -6/ fc])^, (85) 

and in similar fashions for £g and Then, the 4x4 
transformations can be written simply 



bT = 



ilia _ n-l^otvPfvfi 



pa 



0/j 



9 



"/3 



g* B i«*v(i ; (86) 
^^°W^, (87) 



with the 4x4 thermal Bogoliubov and BdG inverse ma- 
trices 



fj,av(3 



-l,a/3 



W; 



a/3 



(89) 
(90) 
(91) 



The existence of the condensate brings no crucial alter- 
ation in defining the thermal counter term, because the 
thermal Bogoliubov transformations Eqs. (f8"0")) and (fSTj) 
remain unchanged. Thus, the thermal counter term is 
the same as Eq. (fTTf 

(\ v-v 
o o ) im ' m 

Note however that the number distribution is that of the 
quasi particles. With the quartet notation, Q can be 
written as 



<5 = 4 £*<(*) «rw 



v 



1 


\ 














1 o) 



The unperturbed and full propagators are given by 

A^O^x') = -t(Q\T[^ a (x)<p^(x')}\Q) , (94) 
G^(x,x') = -i(0\T[^ a (x)^(x')}\0) , (95) 
d^(t,t') = -i(0\T[^ a (t)if(t')}\0) , (96) 

tir p (t,t) = -lim&yitvw) . (97) 



Following the renormalization condition for the non- 
condensed system in Eq. (f3"2")> . which has successfully led 
to the transport equation, we propose the renormaliza- 
tion condition for the condensed system as 



g\j 21 {t 1 t)=Q. 
Let us consider the Dyson equation 



9 11 9 12 



(98) 



y 1 9 22 





Id 11 


° 


)- 


V o 


d 22 ) 




s 11 s 12 
s 21 s 22 



V 1 g 12 
y 1 g 22 , 



, (99) 



where the superscripts denote the BdG indices, a and (3. 
Every matrix elements in the above equation are 2x2 
matrices with the thermal indices, /i and v, which are 
implicit for conciseness of notation. Solving the Dyson 
equation for g , we obtain 



with 



S 



li 



g 11 = d 11 + d^Sg 11 



S 12 (l _ d 22 s 22yl d 22 s 21 



(100) 



(101) 



Obviously, S and g 11 in Eq. (|100[) are upper triangular 
matrices, corresponding to the structures of g^ v and S^ v 
for the non-condensed system. Then from Eq. (jlOOp fol- 
lows 

S# 2 UO=£ fdsds' gl%\t,s)S% m ,{s,s>)g^{s>,t') . 



(102) 

Similarly as in the non-condensed case Eq. (|33[) . let us 
consider the leading order. We approximate the full prop- 
agators g^ 1 ^ 1 by the unperturbed ones eP 1 ^ 1 and 5 12 by 
S 1121 in the right hand side of Eq. (|102|) . Since the con- 
tribution of the thermal counter term to the self-energy 
becomes 



S%%f(t,lf) = -iht(t)6{t-1?)5u> 



{0 1 




V 



\ 



o 

1 oy 



we obtain the following transport equation from Eq. 



(103) 



n € (t) = 2Re / dsS^ op (t,s)e 



iuj£ (t — s) 



(104) 



The tensor form introduced in the previous section 
makes the perturbative calculation very simple and sys- 
tematic. The unperturbed propagators are written as 
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-i6(t-t') 
+ i9(t' -t) 







■ a , s.0 
1 1 



1 1 o 



-iutit — t') ■ 



-iu)t{t — t ) _|_ 



.>n°iw 



i 

( 



6(t - 1') 



V 



+0(t' - 1) 



ll \l + n e {t') 
if \ n e (t>) 

lV( n e (t>) 
1 \l + n e (t') 



l + n<(t)[ |1) 



l + ni(t)\ 111 



y^C 33 )? \vt{ x ') 



i\ l p -iu e {t-t') 







z t {x)\ \zt{x')\ e 



ye(x)\ ly%{x')\ e 



itug(t—t') 



The one- loop self-energy, indicated in Fig. [2] (a), is given in the tensor form as 



9{t-t') 
+0(t' - 1) 
+9(t - 1') 
+9(t' - 1) 
+9(t - 1') 
+9(t' - 1) 
+9(t - 1') 



z e (x)} \z* t {x') \ e 



, 

>\ I iue(t-t') 



\\ rtT'i m {^} Kt -° 



n ei (l + n e2 ) I ll 
(l+n^Jn^J [1 



(1 + n,£ 1 )m 2 I 111 
n ei (l + ne 2 )j [lj 

ll ( U^n^ 

if 1(1 + ^X1 + ^) 



Xzy {x) \ \xly{x') \ e^-^-t') 



(105) 



(106) 



(107) 



min(t,t') 
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FIG. 2: (a) One-loop self-energy diagram, (b) A part of the one-loop self-energy diagrams for a = f} = 1. There are nine terms in total, 
(c) A description in the tensor form for a = j3 = 1. The diagrams can be factorize in terms of tensor form. 



where 



*yy 


= Cvlvl - 




H a vlyf 2 : 


(108) 


Ayz 


= Cvtzt - 






(109) 


Y 

Azy 


= C*tvt - 






(110) 


A.ZZ 


/-a a a 




K Q ^4> 


(111) 



with C, a (x) 



U*(*), 



and a denotes a = 2,1 for 



a = 1,2, respectively. Since 

x Ef»f (x,x')W I - 1 ''"(^)B I " 1 '' 2 (f), (H2) 
we obtain from Eq. (|104|) 



n<(f) = 4 3 2 \ 



di' Re 



e i ^-^i- w «)(*- t ')|( W) Xyy )f {(1 + n/Jn^n^ - n e (l + n h )(l + n t2 )} 

+ e i ^-^+^ t - t \y hXyz )\ 2 {(1 + n t )n u {\ +n ta )- n e (l + n £l )n e2 } 

+e i{ue+ ^- u ' e ^ t - t ' ) \(ye,Xzy)\ 2 {(1 + n t )(l + n ei )n t2 - n e n h (l + m 2 )} 
+e i( w< +a» <1 + W4a )(t-t')|( W) ^)|2 {(1 + nt){l + nei){l + nt j „ n ^ n< j 



r 



(113) 



with 



{ye,x) 



d 3 x y* t ' 



(x) X a (x) 



(114) 



This is the non-Markofhan transport equation for the 
condensed system at one-loop level. The first term in 
Eq. (|113p corresponds to the Beliaev damping and its 
inverse process, and the second and third terms do to 
the Landau damping and their inverse processes. These 
processes bring the system to the equilibrium. On the 
other hand, the fourth term corresponds to a process 
in which three quasi particles are created or annihi- 
lated. We call it the triple production process. The 
existence of the process prevents the system from equi- 
librating, because the collision term is always nonzero 
(1 + rii)(l + u^Xl + ng 2 ) — ngni 1 ng 2 > . Note how- 
ever that if all the energies of quasi particles are positive, 
the process is forbidden due to the energy conservation. 



Conversely, once a negative energy mode exists, the ex- 
citations to the mode will proceed until the condensate 
decays. This exactly corresponds to the scenario of the 
Landau instability in terms of the kinetics. Thus the 
triple production term is interpreted to induce the decay 
processes leading to the Landau instability. 

To see the importance of quantum field theoretical 
treatment based on the proper quasi particle representa- 
tion, let us treat the same system in the particle picture 
of original atoms. Namely, we suppose that a macro- 
scopic number of atoms is in the particular state I = £ c , 
and consider it in the transport equation Eq. (|50|) derived 
for the non-condensed system. The energy ui£ is shifted 
in such a manner that ui a is vanishing, and the particle 
distribution ng(t) is replaced by ne(t) + Se c iN c . If it is 
assumed that N c ~ TV with the total atom number N, 
the transport equation approximately becomes 
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lit 



(t) ~ Ag 2 N c / dt' Re 



+e i( U( -^ 1+Ua )(i- t ') C4t . w {(1 + n< ) n<i (i + n £2 ) - n< (l + n h )n l2 } 



(115) 



This equation corresponds to Eq. (|113[) . but the triple 
production term is absent. Its absence comes from the 
inadequate choice of particle picture. 

The difference between Eqs. (|113[) and (|115[) is deci- 
sive when the system has the Landau instability. We 
conclude that the Landau instability should be described 
by Eq. (|113[) with the triple production term, because it 
is based on the appropriate quasi particle picture. The 
omission of the triple production term, when the en- 
ergy conservation allows it, would violate the unitarity 
in quantum theory. 



V. SUMMARY AND DISCUSSIONS 

In this paper, nonequilibrium TFD is applied to the 
systems of trapped Bose atoms, and the quantum trans- 
port equations of non-Markoffian type have been derived 
for both the non-condensed and condensed systems. We 
have diagonalized the unperturbed Hamiltonians, each 
of which corresponds to the quasi particle picture, and 
this diagonalization procedure in the interaction picture 
is essential for TFD as well as for ordinary quantum field 
theory. To derive the transport equations for the trapped 
systems both without and with a condensate, we have 
applied the self-consistent renormalization condition pro- 
posed by Chu and Umezawa for a homogeneous system. 
In order to make complicated calculations of the self- 
energies transparent, we have also refined the diagram- 
matic calculations in the tensor form, and have developed 
the convenient 4 x 4-matrix formulation in the condensed 
case. Although the transport equations are derived only 
in the lowest order in this paper, their higher order cor- 
rections can be obtained systematically in our method, 
simply by calculating higher order diagrams. In contrast, 
higher order corrections of the transport equations can- 
not be obtained straightforwardly in the other methods. 

For the non-condensed system, the non-Markoffian 
transport equation at two-loop level derived in this pa- 
per becomes very similar in the Markofhan limit to those 
derived in the different methods. While the equations in 
the other methods involve a delta function and require 
a strict energy conservation in each collision, the energy 
in our equation is conserved with the finite width which 
reflects thermal changes and can be calculated. With the 
strict conservation, the collision integral is either zero or 



infinite because of the delta function in a trapped sys- 
tem where the energy spectrum is discrete. Therefore, 
an additional cure was needed to avoid the problem in 
the other methods. Although the problem does not oc- 
cur in the semi-classical method since the delta function 
is integrated over continuous energy variable, the semi- 
classical treatments are not consistent with the particle 
picture in the trapped system. It is remarked that our 
equation in nonequilibrium TFD follows from the correct 
particle picture and needs no additional patch. 

Perturbative calculations are much more intricate for 
the condensed system than for the non-condensed one: 
In the former the two component eigenfunctions of the 
BdG equations complicate expressions. We have merged 
the thermal doublet and the BdG one into a quartet and 
have constructed the 4 x 4-matrix formalism, which is 
helpful for our concrete calculations of nonequilibrium 
TFD. 

In principle similarly as in the non-condensed case, we 
have derived the quantum transport equation in the con- 
densed case at one-loop level. A crucial point in our 
equation is that it involves an additional collision term, 
absent in the transport equations of the other methods. 
Usual collision terms in the lowest order are only those 
corresponding to the Beliaev and Landau damping and 
their inverse processes. Our additional term represents 
creation or annihilation of three quasi particles, and pre- 
vents the system from equilibrating if a negative energy 
mode exists and is suppressed otherwise. So the behavior 
of the system with the Landau instability, described by 
our equation, is distinguished from those under the equa- 
tions without the triple production term. We emphasize 
that the additional term disappears in the inadequate 
particle picture, as was shown at the end of Sec. IV, but 
that it appears naturally in quantum field theory. 

As for the Landau instability, the authors of Ref. (3f| 
have pointed out that the sign of the Landau damping 
rate changes to minus when the system has Landau in- 
stability, which is interpreted as an indication of the de- 
cay of the condensate. In their analysis, the nonequi- 
librium distribution function was roughly approximated 
to the Bose-Einstein one. However, this approximation 
is invalid for the negative energy spectrum, because the 
distribution function becomes negative and therefore un- 
physical. In our scenario with the transport equation, 
the dominant term in the Landau instability is obviously 



11 



the triple production one which is always positive. 

We comment on a relation between the closed time 
path formalism (CTP) [l| and our nonequilibrium TFD. 
Both derive the very similar Dyson-Schwinger equations, 
indeed the same in form. This is because they follow from 
the common Heisenberg equations and Feynman dia- 
grams. It is important to notice that the the same Dyson- 
Schwinger equation does not always give the same solu- 
tion. A main difference between the two approaches is in 
their unperturbed propagators. The unperturbed prop- 
agator in CTP is evaluated over the density matrix at 
the initial time to, out of equilibrium, and would become 

AgJpfoaO ~ [ne(t )]dfy' (t, t')B v ' v [m> (to)] in- 

stead of Eq. (HD in TFD. Thus, while the thermal Bo- 
goliubov matrix in the unperturbed propagator of CTP 
is time-independent and carries information only on the 
initial thermal state, our propagator with the time- 
dependent Bogoliubov matrix adopts temporal thermal 
changes of the system. So the unperturbed representa- 
tion in TFD contains some non-perturbative effects in 
perturbative calculations of CTP. The time-dependence 
of the thermal Bogoliubov matrix in TFD is impor- 
tant also in respect of the self-consistent renormalization, 
lacking in CTP so far: The time-dependent thermal Bo- 
goliubov matrix creates the thermal counter term Q in 
the interaction Hamiltonian, and the derivation of trans- 
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port equations would be impossible without it. 

The transport equation derived in this paper can de- 
scribe only the initial stage of the condensate decay with 
the Landau instability. It is out of our present formula- 
tion to describe a full time evolution of the decay as well 
as that of quantum phase transition through evaporative 
cooling, because we have ignored the time dependence 
of the condensate involving the time dependent quasi 
particle representation. In such cases, the zero modes 
of the BdG equations play a crucial role and must be 
considered. The description of the decay with the Lan- 
dau instability, that with the dynamical instability a full 
description of quantum phase transition are challenging 
subjects in nonequilibrium TFD. 
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